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SUMMARY 

By  an  extension  of  Eward*s  '’diaphragm"  concept  outside  the  wing 
tip,  the  present  paper  presents  two  approximate  methods  for  calculating 
the  aerodynamic  behavior  of  harmonically  oscillating,  sweptback  finite 
wings  with  both  supersonic  leading  and  trailing  edges.  Both  methods 
are  based  upon  the  fact  that  the  contribution  of  the  so-called  "diaphragm" 
on  the  potential  at  a  fixed  point  on  the  wing  is  approximately  canceled 
by  the  contribution  of  a  portion  of  the  wing  itself,  if  the  wing  under¬ 
goes  steady  harmonic  oscillation  at  a  low  value  of  the  frequency  param- 

eter  — -  (where  cd  is  angular  velocity  per  second;  c  is  chord;  a 
a32  - 

is  velocity  of  sound;  and  3  =  \/M^  -  1,  M  being  Mach  number).  It  is 
further  necessary  that  the  point  be  influenced  by  the  diaphragm  from 
only  one  wing  tip. 

The  first  method  expresses  the  integration  in  terms  of  Fresnal 
integrals.  The  final  expressions  for  the  pressure  coefficient  are  in 
the  form  of  integrals  which  may  be  readily  evaluated,  using  numerical 
integration  with  available  mathematical  tables.  Some  calculations  are 
shown,  and  the  results  check  well  with  the  results  of  the  second  method. 

It  is  expected  that  this  method  should  give  accurate  results  in  the 

(DC 

engineering  sense  if  the  frequency  parameter  — x  is  of  the  order  one- 

a32 

half  or  less.  Even  if  the  frequency  parameter  exceeds  this  value,  the 
pressure  increments  calculated  asymptotically  approach  the  correct  values 
at  the  inner  edges  of  the  tip  Mach  cone. 


As  a  byproduct  of  the  above  approach,  the  expressions  for  the  down- 
wash  of  the  wing  wake,  influenced  by  only  the  two  tip  diaphragms,  are 
given  in  terms  of  Fresnal  integrals. 

The  second  method  is  similar  to  the  first  one  except  that  a  series 
expansion  of  the  source- strength  function  is  used.  When  the  frequency 

/  q\ 

parameter  is  small  [ — -  <  —  ,  the  few  terms  used  give  good  engineering 

\a32 
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accuracy.  In  the  present  case,  five  terms  are  considered.  Of  course, 
the  number  of  terms  must  increase  as  the  frequency  parameter  increases. 
This  method  gives  closed  expressions  for  the  pressure  coefficients.  In 
applying  the  method,  the  calculations  are  lengthy,  but  no  great  computing 
skill  is  required.  For  the  benefit  of  design  engineers,  some  expressions 
of  the  average  pressure  coefficient  over  the  entire  rectangular  wing  tip 
influenced  by  three-dimensional  flow  are  given  in  a  closed  form,  so  that 
they  can  be  used  to  estimate  the  tip  effect  in  preliminary  design. 

A  number  of  graphs  are  presented  to  show  the  pressure  distribution 
of  the  wings  at  different  Mach  numbers,  sweepback  angles,  and  frequency 
parameters.  Comparisons  of  the  pressure  coefficients  calculated  from 
the  present  methods  and  the  exact  two-dimensional  solutions  based  on 
linearized  theory  are  made. 


INTRODUCTION 


The  aerodynamic  behavior  of  a  harmonically  oscillating  finite  wing 
in  a  supersonic  flow  has  been  recently  studied  by  a  number  of  investi¬ 
gators,  particularly  Garrick  and  Rubinow  (reference  l),  Eward  (refer¬ 
ence  2),  and  Miles  (reference  3)-  From  the  viewpoint  of  mathematical 
physics,  it  is  one  phase  of  the  mechanics  of  wave  motion  investigated 
by  a  number  of  great  physicists  and  applied  mathematicians  such  as 
Huyghens  (reference  4),  Maxwell,  Kirchhoff,  Sommerfeld,  Hadamard, 
Volterra,  and  others  over  a  number  of  years.  The  classical  treatment, 
however,  is  principally  along  the  line  of  Cauchy's  initial-value  problem, 
while  the  present  aerodynamic  problem  is  a  boundary- value  problem.  The 
classical  results  are  of  some  help  in  the  formulation  of  the  present 
problem  but  not  of  much  help  in  the  solution. 

With  the  advent  of  guided  missiles  and  aircraft  operating  in  the 
supersonic  range,  trouble  has  been  experienced  with  flutter  of  the  wing 
or  control  systems.  It  was  such  troubles  that  spurred  the  author  to 
obtain  an  understanding  of  the  aerodynamic  behavior  of  a  harmonically 
oscillating  finite  wing  at  supersonic  speed.  The  two-dimensional  case 
of  an  arbitrary  unsteady  motion  has  been  solved  in  the  last  few  years 
(references  5,  6,  7,  and  8).  But  in  the  three-dimensional  case  up  to 
now  very  few  useful  results  have  been  achieved. 

The  present  paper  treats  the  case  of  a  finite  wing  with  the  leading 
edge  ahead  of  the  Mach  line.  The  wing  may  be  swept  back  or  forward. 
Another  condition  has  to  be  imposed,  that  the  Mach  cones  at  the  nose  and 
the  tip  never  intersect  within  the  wing  area.  With  slight  modification, 
the  method  can  be  extended  to  more  general  wing  plan  forms.  The  wing 
is  supposed  to  undergo  harmonic  oscillations  of  small  amplitude  during 
flight  at  supersonic  speed. 
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For  the  wing  plan  form  as  specified  above,  the  three-dimensional 
flow  fields  are  limited  to  three  portions  of  the  wing  area:  the  area 
within  the  nose  Mach  cone,  and  the  two  areas  within  the  tip  Mach  cones. 
In  the  remaining  wing  areas,  the  aerodynamic  behavior  can  be  calculated 
with  the  two-dimensional  theory  of  swept  wings  with  supersonic  leading 
and  trailing  edges  under  harmonic  oscillation  (reference  6,  e.g.). 
Therefore,  the  present  investigation  is  limited  to  the  tip  area  and  nose 
area  of  the  wing. 

The  author  should  express  his  appreciation  to  Miss  V.  O'Brien, 

Miss  P.  darken,  and  Mr.  B.  T.  Chu  for  their  assistance  in  carrying  out 
the  project.  The  present  project  was  conducted  at  The  Johns  Hopkins 
University  under  the  sponsorship  and  with  the  financial  assistance  of 
the  National  Advisory  Committee  for  Aeronautics. 


LINEARIZED  DIFFERENTIAL  EQUATION  AND 
SOLUTIONS  OF  SOURCE  INTEGRATION 
General  Analysis 


* 


The  linearized  differential  equation  of  unsteady  compressible  flow 
is  well-known  (references  1  and  5) 


25*  +  »S4  +  . -i/i 'jL  +  u|_ft 

dx^  Sy^  a^\^  ^x/ 


(la) 


or,  briefly. 


(lb) 


where  U  is  the  free-stream  velocity  and  a  is  the  sound  speed  of  the 
free-stream  condition.  (Symbols  used  in_this  paper  are  defined  in 
appendix  A.)  The  velocity  vector  q  =  Ui  +  qf  =  Ui  +  grad  0  and  V2 
is  the  Laplacian  operator.  The  elementary  solution  of  the  potential  at 
a  point  (x,y,z)  at  time  t  due  to  a  source  of  time -dependent  strength  A 
is 

* 


0O(  w,t;  e,n,£>  -  Z  Ti)  +  -  Tg) 


(2a) 
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where  the  source  is  located  at  (!,n, £)  and 


r  =  \/(x  -  i)2  -  32[(y  -  n)2  +  (z  -  (3) 


T  =  (x  ~  g)M  _  L 

1  ap2  .  a 

T  =  iiLlJ l)M  +  £ 

2  a|32  a 


(4) 


where 


In  the  case  of  harmonically  oscillating  source  strength  the  source- 
strength  function  can  he  written  as  a  product  of  two  functions,  one 
containing  the  source  position  and  the  other  showing  the  time  dependence, 
or,  explicitly, 


AU,T),e;,t  -  t)  =  a( s,T|,(;)f(t  -  t) 


Equation  (2a)  then  becomes: 


0o(x,y,z,t,-  i,n,S)  =  -A— jf(t  -  tx)  +  fjt  -  T2jj  (2b) 

In  the  whole  domain  surrounding  the  wing  in  supersonic  flight,  there 
is  one  particular  streamline  surface  which  is  of  interest  to  the  problem. 
This  streamline  surface  is  formed  by  the  streamlines,  which  coming  from 
upstream  either  strike  the  wing  leading  edge  or  pass  through  the  fore¬ 
most  Mach  line  of  the  tip  Mach  cones.  It  follows  the  contour  of  the 
wing  (both  upper  and  lower  surfaces),  reunites  at  the  trailing  edge,  and 
then  extends  to  infinity  at  downstream. 
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Figure  1  shows  the  streamline  surface  which  is  associated  with  an 
oscillating  wing  in  a  supersonic  flow  at  time  t. 

The  whole  streamline  surface  can  he  divided  into  eight  zones  as 
shown.  Ahead  of  the  supersonic  leading  edge  and  outside  the  two  tip 
Mach  cones,  the  flow  is  not  influenced  by  any  of  the  disturbances  caused 
by  the  wing.  This  is  called  zone  1  as  shown.  The  portion  of  the  stream¬ 
line  surface  behind  the  supersonic  leading  edge  and  between  the  tip  Mach 
cones  is  called  zone  2.  In  this  zone  the  potential  at  any  point  on  the 
top  wing  surface  (or  bottom)  is  never  influenced  by  the  flow  on  the 
opposite  side.  Therefore,  the  source  solution  of  equation  (2b)  can  be 
applied  without  restriction,  if  the  source-strength  function  A  is 
known. 

Zone  3  represents  the  portion  of  the  wing  influenced  by  the  right 
tip  Mach  cone  and  zone  4  by  the  left  tip  Mach  cone.  In  these  two  zones, 
the  flow  on  one  wing  surface  is  influenced  by  the  flow  on  the  opposite 
side.  The  simple  source  solution  of  equation  (2b)  cannot  be  applied. 

The  present  paper  gives  the  approximate  solution  for  these  two  zones 
for  the  case  where  the  two  tip  Mach  cones  meet  on  or  outside  the  trailing 
edge.  The  present  method  with  slight  modification  can  also  be  extended 
to  the  case  that  the  tip  Mach  cones  intersect  inside  the  wing  area  (not 
shown  in  fig.  l),  but  no  diaphragm  within  the  forward  Mach  cone  of  the 
point  under  consideration  is  influenced  by  the  diaphragm  on  the  other 
side. 

Zone  5  is  the  portion  of  the  streamline  surface  outside  the  wing 
surface  influenced  only  by  the  right  tip  Mach  cone;  similarly  zone  6  is 
influenced  only  by  the  left  tip  Mach  cone.  Zone  7  is  the  portion  of  the 
streamline  surface  after  the  trailing  edge  influenced  by  neither  Mach 
cone.  Zone  8  is  the  portion  of  the  streamline  surface  influenced  by  both 
Mach  cones.  Of  course,  if  the  wing  moves  slightly  up  and  down  with  time, 
the  streamline  surface  and  consequently  the  diaphragm  must  also  move  up 
and  down  with  it.  With  the  concept  of  Eward  (reference'  9),  any  finite 
wing  in  supersonic  flow  can  be  considered  to  be  attached  over  its  stream¬ 
line  surface  outside  the  wing  (such  as  zones  5>  6,  7,  and  8)  to  a  thin 
weightless  rubber  diaphragm,  of  which  the  pressure  on  both  sides  at  a 
point  must  be  balanced  and  through  which  no  fluid  can  leak. 

The  diaphragm  is  imagined  to  extend  to  the  outer  edge  of  the  tip 
Mach  cone,  and  to  coincide  with  the  shape  of  the  streamline  surface. 

Thus,  the  rubber  diaphragm  can  be  considered  as  an  extension  of  the  wing. 
The  flow  on  the  upper  side  of  the  wing  with  its  extension  can  be  con¬ 
sidered  as  independent  of  the  flow  on  the  lower  side  of  the  wing  and 
vice  versa.  The  potential  at  any  point  on  one  side  of  the  wing  can  be 
evaluated  from  the  source  distribution  on  the  same  side  of  the  wing  and 
the  diaphragm  within  the  forward  Mach  cone  of  that  point.  In  the  steady 
case  it  has  been  shown  by  Puckett  that  the  source  strength  at  a  point 
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n 

of  the  wing  is  directly  proportional  to  the  slope  at  the  point  and  inde¬ 
pendent  of  the  slope  of  its  neighboring  points.  Eward  has  shown  that  * 

the  solutions  obtained  using  the  diaphragm  concept  are  identical  to 
those  of  conical  flow  (reference  9) • 

For  steady  flows,  the  principal  advantage  of  Eward '  s  method  is 
that  the  slope  of  a  point  on  the  diaphragm  is  determined  solely  by  the 
slopes  of  the  wing  surface  along  the  characteristic  or  Mach  line  through 
the  point,  as  long  as  the  point  on  the  diaphragm  is  ahead  of  the  Mach 
cone  from  the  other  wing  tip.  In  the  case  of  unsteady  flow,  it  can  be 
shown  that  the  above  statement  is  approximately  true  if  the  wing  slopes 
vary  very  slowly  with  time.  Using  the  diaphragm  concept,  the  disturbance 
potential  at  any  point  can  be  considered  as  the  effect  of  sources  of  the 
wing  system  on  the  same  side.  (Of  course,  the  slope  of  the  imaginary 
diaphragm  depends  on  the  flows  over  both  wing  surfaces.)  It  becomes 
apparent  that  the  potential  can  be  evaluated  by  the  well-known  source- 
integration  method.  Let  us  consider  a  point  on  the  diaphragm  where  the 
potential  contributed  by  sources  on  either  side  of  the  wing  must  be  the 
same.  It  has  been  shown  by  Eward  that  the  source  strength  is  propor¬ 
tional  to  the  vertical  component  of  the  perturbation  velocity,  or,  more 
directly,  is  proportional  to  the  local  slope  of  the  wing  surface  or  the 
diaphragm  even  if  the  flow  is  not  steady.  Mathematically,  this  is  shown 
«  w  Ua 

as  A  =  ^  =  5  where  w  is  the  vertical  component  of  the  pertur-  m 

bation  velocity  and  a  denotes  the  slope  of  the  streamline  surface. 

* 

The  following  equations  can  be  written  to  denote: 

The  slope  at  the  upper  wing  surface  at  time  t  -  by 

aT,l  =  °T(^Tl>+0>'t  "  Ti) 

The  slope  of  the  upper  wing  surface  at  time  t  -  Tg  by 

°T,2  -  -  t2) 

The  slope  of  the  lower  wing  surface  at  time  t  -  by 

-  ti) 

The  slope  of  the  lower  wing  surface  at  time  t  -  by 

°B, 2  =  ~  t2) 

The  slope  on  the  diaphragm  at  time  t  -  by 

h  =  H6,T),0,t  ‘  Tl) 

The  slope  on  the  diaphragm  at  time  t  -  by 

X2  =  x(e,ri,o,t  -  t2) 
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Then  the  potential  can  be  written  at  any  point  on  the  top  or  bottom 
streamline  surface  as: 


0T(x,y,+°,t) 


U  P  P(a T,1  +  ctT,2)  d^  dT> 

^ lJJ  7 


w 


^  Xp  +  ^2^  d£  dr) 
r 


'’N 


"°^t) 


(6) 


u 


Or,  in  an  oblique  coordinate  system  whose  axes  lie  parallel  to  the  Mach 
lines  and  whose  origin  is  placed  at  the  point  of  tangency  of  the  Mach 
line  to  the  leading  edge  of  the  wing  tip  (the  junction  of  the  supersonic 
and  subsonic  leading  edge)  (see  fig.  2),  these  equations  can  be  written 
as: 
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0T(u,v,4O,t)  = 


aT(u',v',t 

'  “  Tl) 

I  +  CTrpI 

(u',v',t  -  t2) 

V 

((u  -  U '  )  ( V  ■ 

-  V) 

du1  dv’  - 


u  I  I  x(u’,v',t  -  TX)  +  x(u',v',t  -  t2)  ^  ^ 
2itM 


\j( u  -  u'  )(v  -  v'  I 


0B(u,v,-O,t)  =- 


U 

2m 


Swl+Sw2 


°b(^W  -  Ti)  +  ^(uSvyb 

]f( u  -  U')(v  -  v') 


du'  dv' 


-2-  r/>x(tt'>v'<t  -  t0  tx(u',v,t  -  t2)-'..  av, 

/(u  -  u')(v  -  v') 


2nM 


’D 


(7) 


where  (u,v)  is  the  transformed  coordinate  of  the  point  originally  at 
(x,y)  and  (u',v')  is  the  transformed  coordinate  of  the  source  element 
originally  at  (|,t|)  in  the  x,y  system.  The  old  and  new  coordinate 
systems  are  related  by  the  following  transformation  equations: 


% 


w 


w 


w 


u  =  ik  (x " Py) 

X  =  |  (v  +  u) 

u'  =  W  ^ 

6'  =  |  (▼'  +  u') 


v  =  (x  +  Py) 


y  =  |  (v  -  u) 


▼'  =  ^  (6  +  Pn) 


n  =  |  (▼'  -  u') 


(8a) 


d|  dT]  =  ^j|  du'  'dv' 


quantity  r  and  the  time  intervals  t^_  and  T2  in  the  new  coordinate  system  are 


■ 
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0) 
> 
•H 
0)  -P 
U 

■P  (D 

_  Ph 
rtf  CQ 
P  0) 
cti  U 


•2(u') 


(A  similar  expression  can  be  obtained  for  Now  it  will  be  shown  that  in  the  case  of  har¬ 

monic  oscillation  the  second  integral  is  approximately  zero  under  certain  conditions. 
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•^Members  of  the  MCA  staff  recently  expressed  the  opinion  that  "Since  t-j_  and  t 2  con¬ 
tain  the  variables  uT,  a  second-order  integral  equation  is  generally  required  rather  than  the 
simplified  first-order  integral  equation  given.  Equation  [1 3)  is  satisfactory  if  t  may  be 
represented  as  varying  linearly  with  time  in  the  interval  (Ti  “  T2)*"  Preserr^  analysis 

the  use  of  equation  (13)  is  justified. 
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1  =  SF(TD  -  v')  T2 


Substituting  equations  ( 11a) ,  (lib),  and  the  above  value  of  the  t's 
into  equation  (13),  there  is  obtained: 


D 


4(u'  ,v'  )e 


ics 


W) 


v2(u' ) 


.  v*:V2 


dv'  = 


fa  - v') 


'Vo(u') 


V^u*) 


-4-^) 

ae  _ 

(VD  -  v')l/2 


dv'  (14) 


Now  if  in  equation  (12)  cos 


IhK  •  u')1/2fa  -  T’)1/2" 


=  1  exactly, 


the  integral  at  the  right-hand  Bide  of  the  equation  will  reduce  to  the 
form: 


Vo(u')  -i(U 


ap 


N-V) 


a 


'vi(u') 


(vw  -  V')1/2 


W  — /v  -v') 

dr.  ♦  |  l 

Wu’)  fa-fal/2 


dv 


(15) 


which  according  to  equation  (l4)  is  zero,  because  here  vw  =  v^ 


In  general,  cos  ^(u^  -  u')1/2^  -  v')1/2 


=  cos 


mr 


is  less 


than  unity.  (Note  that  r  =  >/(x  -  £)2  -  p2(y  -  q)2  is  the  hyper - 

v 


bolic  distance,  not  the  geometrical  distance.)  But  for  small  values  of 
CUT  mr 

— ,  the  variation  of  cos. —  from  unity  is  very  small.  Suppose  that 
a  a 

the  following  condition  is  imposed:  —  <  Then 

05T  -  1  /  mr  \2 


1  >  COS 


■3^1  -7:  =  0.875.  In  both  integrals  the  factor  cos  — 

a  \  a  /  a 

will  reduce  the  value  of  each  integrand  a  small  percentage  in  view  of 

oxr  nr 

—  for  small  values  of  — . 
a  a 

Of  course,  the  net  contribution  of  both  integrals  cannot  cancel  out 

cur 

exactly  as  if  cos  —  =1,  which  is  located  on  the  left  tip  Mach  line. 

a 


the  monotonically  decreasing  nature  of  cos 
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Obviously, 

X  — >  oo  as  Vj-)  — >  v2 
X  — >  0  as  vD  — oo 


For  the  case  vD  >  v2  (fig-  2),  v2  -  Vj  can  be  divided  into  n  equal 
parts  and  the  ith  division  can  be  denoted  by  v^  (where  vi'  = 
for  i  =  0  and  =  v2  for  i  =  n) .  Then, 


(19) 

where  the  division  is  so  fine  that  v2  -  v'  =  v2  -  v^  can  be  taken  out 

of  the  integration  sign  for  the  first  approximation.  Incidentally,  the 
integral  left  can  be  transformed  to  the  Fresnal  integral  (reference  10). 
Equation  (19)  then  becomes : 


C(Zi)  -  iS(Z1+1)  +  iS(Zi) 


(20) 


NACA  TN  2b6rJ 


l6 


where 


Zi  = 


(  2(0  \l/  2 

Up*/ 


c(z) 


£ 

2 


Z2  dZ 


S(Z) 


sin  |  Z2  dZ 


(21) 


Some  of  the  basic  properties  of  the  Fresnal  integral  are  given  in 
equations  (29)- 

APPROXIMATION  THEORY  OF  VELOCITY  POTENTIAL  AND  PRESSURE  COEFFICIENT 
AT  A  POINT  ON  A  HARMONICALLY  OSCILLATING  SWEPTBACK  FINITE 
WING  WITH  SUPERSONIC  LEADING  EDGE3  -  FIRST  METHOD 

Assume  the  whole  wing  flapping  up  and  down  according  to 
a,p(t)  =  aeiayk  where  a  is  the  amplitude  of  oscillation.  In  fact. 


Maximum  descending  rate  of  wing  h^x 
Free -stream  velocity  U 


The  frequency  of  oscillation  per  second  is  — .  The  velocity  potential 
of  such  oscillation  may  be  considered  as  due  to  harmonic  sources  and  is 


^The  author  should  express  his  appreciation  to  Dr.  J.  C.  Eward 
for  a  discussion  on  this  approximation  method  in  May  19^3. 


given  by  equation  (6).  It  has  also  been  shown  that  for  the  harmonically  oscillating  wing 
equation  (6)  can  be  further  simplified  to  equation  (l6)  as  follows: 
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Since  in  most  practical  cases  straight  leading  edges  are  under  consider¬ 
ation,  put  (cf.  equation  (9)  and  fig.  3): 

v  =  V]_(u)  =  k^u  or  u  =  u-^(v)  =  v/k^ 
v  =  v2(u)  =  l^u  or  u  =  u^v)  =  v/k2 


as  the  equation  of  the  supersonic  wing  edge  outside  the  tip  Mach  cone 
and  the  subsonic  leading  edge  within  the  Mach  cone.  Then,  equation  (22a) 
can  he  written  as: 


0(uw>vv>+o>t)  = 


...  liot 
Uoe 


_  du' 

2Mit  /  (u*  -  u')1/2 

JVk2 


v 


kiu 


-Ixo 

dv'  a3 

172  6 


[(vw“v,-fuw‘ul )  +  |(uw-u'  )^2(  vw‘v’  )1//2] 


(Yw  -  V') 


(22b) 


where  it  is  understood  that  e 

stands  for  the  sum  of  e 


-io) 

ap 


[(vw-v'+uw-u' )  +|(uw-u'  )1/2(vv-v’  )ly/2] 

^[(vw-v'+uw-u' )  -  §(uw_u'  )1/2(vw-v'  )ly/2] 

Tjf  (( vv- v  ’  +uw"u ' )  + f(  uw"u ' ) 1/2  ( vw- v ' )  1/2J 

and  e  ~J 

Introducing  the  following  new  variables  Z  and  Z  instead  of  v' , 

2  ■  (J£)1/2[(t''  -  T')l/2  -  5  (“»  -  "')1/2]  <«“) 


2  ’  (^)1/,2[(vw  -  T')1/2  +  j?  "  u')1/2 


(23b) 


bo  that 


=  dz  = 

*)  Sir  _  v'N1/2 


dZ  =  dZ 


2  \a0 


(vv  "  v') 


Then  equation  (22a)  can  he  written  as: 
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Nov  introduce  another  new  variable  Y  instead  of  u' : 


so  that 


1  /2off  \ V 2  du' 

2 VitaM2/  (uw  -  u')1/2 


Z 


1  ~ 


2m 

apn 


kluv) 


Y 

P 


(26) 


(27a) 


(27b) 


2co 

a0Jt 


(vv 


*1"*) 


(27c) 


(27d) 


If 


Y1  = 


_  /aoB\i/2/  _  M1/2 


naM" 


^uv 


k, 
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where  j(  )  is  the  Bessel  function  and  )  is  the  Hankel  function 

of  the  second  kind  and  various  subscripts  indicate  the  order  of  these 


functions.  For  large  values  of  e, 

c<«>  ■  I +  h sln  I  +  °(^)  <29f> 

S(€)  =i-icoS|e2  +  0(i)  (29g) 

» 

S(±°°)  =  C(±~)  =  ±°°  (29h) 

S(e)  =  -S(-e),  C(e)  =  -C(-e)  (29i) 


For  details,  see  reference  10,  pages  5^  to  and  7^  to  7^5  >  and  also 
reference  11. 

However,  it  is  difficult,  if  not  impossible,  to  reduce  equation  (28) 
to  a  closed  form.  Besides,  there  is  usually  more  interest  in  the  pres¬ 
sure  coefficient  which  is  generally  known  as: 


Cp(x,y,t) 


_2_^  _  2  ^ 
U2  St  U  Sx 


(30) 


rather  than  the  potential  itself.  In  the  u,v  coordinate  system,  the 
pressure  coefficient  is  given  hy: 


Cp(x,y,+0,t) 


_2_  ^  _  2/S0  duw  +  Svv\ 
U2  St  U\Suw  Sx  Svv  Sx/ 


_2_  ^0  M  +  S0  \ 

U2  St  UPySuw  +  Svw// 


(3D 


Now, 


80 

8t 


(32) 
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Differentiating  equation  (28)  with  respect  to  u^.  nnH  .  vw. 


_i£  y  2 

UScle^  e  2  1 


8u 


■w  itM2 


-C(Z3)  -  C(z3)  +  i[s(z3)  +  S(z3| 


^  + 


PY1  .iiY2  /-|zi2  -fzi2' 

e  2  a\£-l -  +  O— 

\  PZX  +  Y  PZ-L  -  Y 


(33> 


80  =  __1_  upae^  e _ f_ 

8vw  kg  nM2  Y^ 


-^Yx2 


[*(Z3)  -  c(z3)  +  1S(Z3)  +  1S(Z3)]  - 


Uoe1^ 


P*1  i* 


-  4r  Y2 


_  ill  7  2  7  2 

2  *  |  "  2  Z1  -TZ1 

e  dYl- -  +  - 

PZX  +  Y  PZX  -  Y 


(34) 


where 


Z3  ■" 


2o)  y  ,  k]M^  p 

apjt  (Vw  "  lUw)  +  ~ p2~  Y1 


1/2 


Yl 

P 


(35a) 


z3  = 


2cjd  kiM  / 

—  (VW  -  *l«v)  -  -p-  Yx 


1/2 


(35b) 
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The  pressure  coefficient  is  then  given  by: 


Cp(x,y,*),t)  =  e  2  ^  [-C(Z3)  -  C(Z3)  +  1S(Z3)  + 


(z3)  + 

J  «Jo 


1  _  1*  y2  f 

dYe  2  rf  Fc(zi)  ■  C(^) +  iS(zi) + 


(1  -  kl} 


i*  2  i£  7  2 
-  ~  Z1  -  ~  Z1 

- - - —  +  £_ - 

,  (3Z1  +  Y  3Z1  -  Y 


The  above  equation  is  used  to  calculate  Cp  along  the  trailing 

edge  of  a  rectangular  wing  tip  which  is  harmonically  oscillating.  (See 
fig.  5.)  In  this  particular  Instance,  M  =  \[2}  the  frequency  parameter 

CDC 

— X  =0.2,  and  the  chord  c  is  taken  as  unity.  The  curve  is  slotted  in 

aP  C., 

the  form  - *L.  against  y  (positive  y  measured  outward  along  the 

-  ioit 
ae 

span)  so  that  for  any  specified  amplitude  of  oscillation  a  the  pres¬ 
sure  distribution  at  any  instant  t  can  be  estimated  from  a  single 
Op 

curve.  As  — jf-r  is  a  complex  quantity,  the  real  and  imaginary  parts, 
a e11013 

denoted  by  (  and  ( — r^r)  >  respectively,  are  plotted  separately. 

\aeiayt7p  \aelayt7j 

If  the  harmonic  oscillation  Is  a  cos  ait,  that  is,  the  real  part  of 
^iait^  then  c  can  be  calculated  from: 


cp(?h  - a  (= 


^^4 


COS  Oit  - 


sin  ait 


(37a) 


Similarly,  if  the  oscillation  is  a  sin  as t,  that  is,  the  imaginary  part  of  ae1^  ,  it  is  given 
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in  equation  (25c) 
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and,  r eplac  ing  k^_  by  k^ 1  , 


Zj"  =  Zj_  in  equation  (25a) 

Zl"  =  Z1  in  equation  (25c) 

Similarly  Y  can  be  defined  as  in  equation  (26)  and  by  the  same 
procedure  C  can  be  determined* 


APPROXIMATE  ANALYTIC  SOLUTIONS  OF  PRESSURE  ON  AN 
OSCILLATING  FINITE  WING  -  SECOND  METHOD 


In  the  preceding  section,  an  approximate  integral  solution  for  the 
oscillating  finite  wing  is  given,  and  an  example  of  calculation  is  shown, 
but  the  full  evaluation  of  the  integral  for  a  large  number  of  cases 
would  take  more  time  than  was  available.  However,  for  very  low  fre¬ 
quencies,  this  problem  can  be  solved  asymptotically  by  using  a  few  terms 
of  a  series  expansion  of  w(u,v,t).  In  the  following  section  the  approx¬ 
imate  analytic  expressions  will  be  shown. 


Analysis 

From  equation  (22a), 


0T(x,y,+o,-t) 


(40) 


where  Sw^  is  the  area  PABC  as  shown  in  figure  .6.  Since  the  two- 

dimensional  solution  with  supersonic  leading  edge  is  known,  the  potential 

can  be  obtained  alternatively  by  calculating  the  contribution  of  S’ 

"  wa  -  * 
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which  is  the  area  BCD,  the  difference  between  the  two-dimensional  source 
area  PAD  and  Sw^  as  shown  in  figure  6.  Then, 


0T(x,y,+°,t)  = 


^-icsTi  +  -ioyr2 
r 


d!  df]  - 


=  0o(x>y^+0^)  -  $(x>y>+°>‘t) 


(41) 


0(x,y,+O,t) 


-iODT  -iO)To 

e  +  e  ^ 


r 


d|  dr] 


1 

np2 

(42) 

The  variables  r-^,  Tg,  r,  and  w  are  the  same  as  before.  In  this 

development,  it  is  found  more  convenient  to  introduce  a  new  coordinate 
system  as  follows:  The  x'-axis  passing  through  the  fixed  point  (x,y) 
runs  upward  along  the  right  forward  Mach  line.  The  origin  is  taken  as 
the  point  where  the  x'-axis  intersects  the  wing-tip  side  edge.  The 
y'-axis  following  the  other  Mach  line  runs  left  upward  as  shown  in 
figure  6.  The  origin  of  the  x',y'  coordinate  system  is  translated 

from  the  origin  of  the  x,y  system  a  distance  where  5  can  be 

shown  to  be 


w(bn^+0,t  -  T-|_)  +  w(|,T],+0,t  -  Tg) 


d|  dT) 


wa 
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5  =  |  (x  +  Py) 


which  is  a  function  of  the  point  P(x,y)  only.  The  transformation 
relations  are 


(43) 


x'  =  - 


M_ 

23 


(x  -  |  5  -  py)  =  -§@f  x  -  My  -  &)  x  =  -|(x’  +  y'  -  &) 


y.  P 

*  2p 


(x  -  £  6  *  Py)  -  -i(|  x  +  mr  -  s)  y  =  i(x'  -  y’) 


n'  =  - 


M 

2P 


(i  - 1 6  +  e,) 


i  =  -|(S’  +  11 '  "  &)f  (44) 


i-h(5’  -i') 


a(x',y')  =  *£ 

d(x,y)  23 


and 


r  =  pm  -  s')  =  ^ vv^7  - x') 


T  -  -L 
1  _  ap 


T2  "  ap 


U’  +  V  -  x')  -  -  yvU'  -  x') 


(!'  +  V  -  x*)  +  -  x') 


(45) 
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Integrating  with  respect  to  q  '  and  introducing  a  group  of  new 
notations, 

PS 

.PC~Trrr_iO)t  PP  r- 


0(x',y',+O,t)  = 


-\/27^Ude' 


^  +  e0z®  + 


Jo  v  +'6Z  +e8_  ’ 

e10zl°)  +  "  z2KzJ|  +  Sgz6  +  g8z8)  +  (b  ■  z2)3//2h6Z6 


where  ^  -  rh>  -  ftp zC  - b  -  f2  *>  *  -  ^  -d 

the  e's,  g's,  and  hg  are  given  in  appendix  B.  With  the  aid  of  the 
formulas  in  appendix  C,  the  following  integration  of  equation  (49)  can 
be  carried  out  with  respect  to  z: 

0U,y,+O,t)  =  -i  ae^  si11'1  l/fZ-ff  H  ffl.n(x  ‘  ty)"  +  \MX  +  3y)  JZ  °o  n(x  ’  + 

v  n=l  n=0  ’ 

fyU  +  ey)3  o3  n(x  -  ky)n  +  ll-y(x  +  py)5  ^  cr,  (x  -  ky)n  + 
n=0  ^  n=0  4,n 

\/-y(x  +  py)7  2Z  a5^ n( x  -  ky)n  +  ]j-y(x  +  py)9  a6^0  + 


2  _  v  _ 


V-y3(x  +  Py)  XI  a7,n(x  -  *y)n  +  \|-y3(x  +  Py)3  £2  n(x  -  fey)n  + 

n=l  n=0  3 

, -  1  _ ____ 

V~y3( x  +  Py)5  ZI  a9,n(x  -  ^y)11  +  l/-y3(x  +  Py)7  a 

n=0  10> 0 


where  the  a^^’s  (m  and  n  are  integers)  in  terms  of  a>,  k,  M, 
and  p  are  given  in  appendix  B.  By  definition, 

C  =  _  2 

P  U2  6t  U  6x 

Differentiating  equation  (50)  with  respect  to  t, 

2  6^  2ico  -r 

U2  dt  U2 


(52) 


Differentiating  equation  (50)  with  respect  to  x  and  adding  equation  (52) 


10,0 
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where  agaih  the  am  's  are  given  in  appendix  B.  The  actual  pressure 
coefficient  is  ’ 

_  2  2  ^T 

P  U2  dt  U  ox 

=  (  2_c#  2 

U2  3t  U  bx  y  ^2  3t  U  bxj 

=  cpo(x,y,0,t)  -  Cp(x,y,0,t)  (54) 

where  0T  =  0O  -  0  as  given  in  equation  (4l)  and  C  (x,y,0,t)  is  the 
pressure  coefficient  of  the  two-dimensional  case  as  given  in  reference  5. 


Rectangular  Wing  Tip 

i 

For  the  rectangular  wing  (k  =  0)  equations  (53)  and  (54)  reduce  to 


1 

,3 

H|  ^ 

sin" 

¥ 

3 

n=0 

+  \/-y(x  + 

2 

V-y(x  + 

Py)3 

Z 

cnxI1  +  \j-y(x  + 

Py)5  Z 

n=0 

1 

n=0 

\f~; y(x  + 

Py)7 

r 

enxI1  +  V-y(x  + 

Py)9  f  + 

n=0 

3 

2 

\|-y3(x  + 

Py) 

z 

n=0 

6nxI1  +  \j-y3(x  H 

h  Py)3 

n=0 

\|-y3(x  +  py)5  +  \/-y3(x  +  py)7 

n=0  ’ 


(55) 


where  the  a's,  b's,  c's,  d's,  e's,  f,  g's,  h'B,  i's,  and  j 
are  given  in  appendix  B.  (The  particular  om  n's  for  k  =  0  are 

found  in  appendix  B.) 
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3^ 


.  .  Average  Pressure  Coefficient  over  Rectangular  Wing  Tip 

For  engineering  purposes,  an  average  value  of  Cp^av  over  the 

entire  tip  influenced  by  the  three-dimensional  flow  may  be  useful. 
(Assume  unit  tip  chord.) 


-p,av 


= 

2P 


■’l  P-x/p 

dx  dyC, 

o  Jo 


i-x/p 


=  2P  '  dX 

10  ^0 


dy 


Jcp(x>y)  -  cp(x  +  0y>+o) 


Cpl  +  Cp2 


(56) 


where  2|3  is'  the  wing-tip  area  influenced  by  three-dimensional  flow. 


'Pi 


-1  -  lent 
2n(3  ae 


dx  ax 
n 

n=0 


n 


p-x/p 


dy  sin" 


^0 


-1  lx  +  (3y  n 

V  x 


.1  =  ol 

1  -  icot  I  ,  —  n 

ae  dx  2_  V 

'o  n=0  J-l 


cos  1  q  dq 


-1  -..icot  <sr-  an 
n=0 


(57) 


where 


x  +  23y 
x 


q 
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CP2  =  dX 


l-x/P  r  4. 

dy  \/-y(x  +  Py)  bnxn  + 


l-y(x~+  Py)3  5^  cnxn  +  J-y(x  +  Py)5  51  dnxI1  + 
n=0  ;  n=0 

_  1  _  _ _  3 

f-y(x  +  Py)?  31  enxI1  +  \/-y(x  +  Py)9  f  +  J-y3(x  +  (By)  31  + 


-y  (x  + 


Py)3  3_  h^x11  +  l/-y3(x  +  Py)5  3_  inxI1  +  \Z-y”3~(x  +”pyj  •  j 
n=0  7  n=0  v 


(58a) 


Let  Py  =  -z^.  Expanding  higher  powers  in  terms  of  (x  - 

the  integral  may  "be  integrated  by  formula  (C2)  in  appendix  C  to 


CP2  -  dx(bi'*2  +  V*3  +  b3'x4  +  V*5  +  »5'x<5) 


1  -  icot  3“  ^n' 

—  cue  >  — — - 


5St  n  +  2 


(58b) 


The  hn'  s  are  given  in  appendix  B.  Hence, 


5 

o  -I  -  lost  sr~  an  +  ^n 

cp,av  ’  JS  “  E  "n+T" 


E  nTT"  *  -=•“  E  bn 

n=0  n=0 


(See  appendix  B.) 


Calculations  and  Graphs 


Owing  to  the  complicated  nature  of  the  formulas  in  this  section^, 
the  calculations  are  rather  tedious.  The  calculations  are  too  long  to 
be  included  here,  but  the ■ essential  results  are  summarized  in  figures  7 
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to  15.  The  pressure  coefficient  given  in  this  paper  is  referred  to  a 
-point  on  the. upper  surface  of  the  wing.  Figures  7  to  10  give  the  pres¬ 
sure  distribution  on  a  rectangular  wing  tip  of  unit  chord  with  different 

free-Btream  flight  Mach  numbers.  The  frequency  parameter  —  is  taken 

*  a 

.  .  •  '  .  , ,t  jtC 

as  0.2  in  these  calculations.  As  before  the  complex  quantity  — 

ae-tot 

is  used  and -plotted  in  two  parts  -  real  and  imaginary  (e.g.,  figs.  7(a) 
and  7(h)).-  These  curves  are  plotted  along  the  spanwise  direction  (y 
being  negative  when. Measured  inward  from  the  tip)  at  various  percentages 

of  chord.  If  —  along  the  chbrdwiee  direction  is  desired,  some 
aeXLUU  ; 

cross -plotting  is  necessary.  It  should  be  noted  that  these  curves  do 
not  give  the  pressure  coefficient  Cp  directly  but  Cp  may  be  calcu¬ 
lated  easily  from  them  by  using  equation  {37a)  or  (37b)  or  both. 


Figure  7  is  calculated  for  a  Mach  number  of  {2 .  (Note  that 
figure  7  checks  well  with  figure  5  which  was  calculated  by  the  first 

method.)  It  is  seen  that  ^  j^t  j  zero  ai  the  wing  tip  (see 

7(^))  for  all  values  of  x,  except  at  the  leading  edge  where 
/  \  '  ; 


[  )  becomes  2it,  which  is  in  agreement  with  the  two-dimensional 

\ae  /r 

result  along  the  Mach  line  -y  =  x.  Taking  the  trailing  edge  (x  =  1.0) 

/  nCL  \ 


as  another  example,  ^  — J  increases  from  zero  at  the  tip  (y  =  0)  to 

a  value  of  -1.885^  at  -y  =  x  =;  1.  Thus  it  is  seen  that  for  this 
/  nCL  \ 


CaSe  J  varies  very  little  along  the  left  Mach  line  of  the  tip 

as  is  shown  by  the  dotted  line  in  figure  7(a).  The  horizontal  solid 
line  for  -y  >  1.0  is  the  two-dimensional  result.  The  imaginary  part 

/  *cp  \ 

l-  last)  shown  in  figure  7(b)  is  of  a  similar  nature.  Figures  8,  9, 

\ae  /j 

and  10  show  the  same  wing  oscillating  under  the  same  conditions  except 
that  the  free-stream  Mach  number  is  now  2,  3>  and  4,  respectively.  A 
few  points  may  be  of  interest.  As  the  Mach  number  increases,  the  maximum 

magnitudes  of  both  ( \  and  \  decrease,  while  the  two- 


aei^ty^R 


dimensional  values  become  more  nearly  constant. 
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In  order  to  show  the  accuracy  of  this  method,  a  comparison  between 
the  result  of  the  present  approximate  method  and  the  two-dimensional 
values  along  the  left  Mach  line  of  the  tip  cone  is  given  In  figure  11. 

The  case  considered  is  a  rectangular  wing  tip  (k  =  0)  and  the  free-stream 
Mach  number  Is  taken  as  )[2 .  The  Solid  line  shown  is  the  approximation 

{JOC 

and  the  dotted  line  is  the  exact  one.  At  — —  =  O.283,  the  deviation  is 

aP^ 

oyc 

unnoticeable.  At  — —  =  O.707,  which  is  beyond  the  imposed  limit*  the 

aP^  •  '  N  .  '  '  '* 

maximum  deviation  of  the  imaginary  part  amounts  to  10  percent.  The  same 
curve  at  x  =  0.8  corresponds  to  the  case  of  x/c  =  1  with  lower  fre- 

QJ  0  X 

quency  parameter  — -  =  0.707>  at  —  -  O.5656  the  maximum  deviation  is 

ap^ 

of  the  order  of  2.  percent.  Thus  beyond  =0.5  the  approximation  of 

ap^ 

the  imaginary  part  is  not  very  favorable  unless  more  terms  of  the  series 
expansion  are  taken. 


There  is  another  parameter  which  influences  the  accuracy  of  the 
present  method,  that  is,  the  sweepback  angle  (A  =  tan~^  k)  of  the  super¬ 
sonic  leading  edge.  Figure  12  offers  such  a  comparison.  Here  the  Mach 


number  is  \j2, 


0.283*  and  the  parameter  k  varies.  The  heavier 


line  gives  the  two-dimensional  exact  results  along  the  left  tip  Mach 

Op 

line,  and  the  lighter  line  of  the  same  kind  gives  — of  the  present 


method .  Although 


ae1^ 


does  show  slight  deviation  from  the  exact  value 


at  higher  values  of  k,  the  comparison  is  favorable. 

In  order  to  show  some  cases  with  high  sweepback  angle  (although  the 
leading  edge  is  still  supersonic),  figure  13  gives  the  distribution  of 
CP 

~"i7ivh  over  the  three-dimensional  wing  tip  of  a  45°  sweepback  angle  at 
ae 

M  =  2.  For  a  matter  of  convenience  in  calculation,  the  results  are 
expressed  -in  oblique  coordinates  x  +  py  and  x  -  y.  The  line 
x  +  Py  =  Constant  is  parallel  to  the  left  Mach  line,  while 
x  -  y  =  Constant  is  parallel  to  the  leading  edge.  Here,  the  chord  * 
normal  to  the  leading  edge  is  chosen  as  unity.  Therefore,  x  -  y 
represents  percentage  of  chord.  With  such  a  definition  of  the  chord. 


the  chord 


used  in  the  frequency  -parameter 


should  be  considered 


as  this  normal  chord.  The  circles  at  the  left  lower  corner  of 
figure  13(a)  represent  the  values  along  the  left  Mach  line.  With  the 
sketch  on  the  figure,  the  general  features  of  the  curves  are  clear. 
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A  comparison  of  results  of  this  sweepback  case  with  two- 

-  '  cB. 

dimensional  exact  values  of  — —  is  also  shown  in  figure  l4.  The 
!  .  ae1(0t 


.  c 

comparison  is  favorable.  The  two-dimensional  value  of  — *— —  for  a 

aeIZD’t 

sweptback  wing  is  obtained  from  reference  6. 


For  engineers,  the  pressure  distribution  on  the  three-dimensional 
wing  tip  is  not  of  so  much  interest  as  the  average  pressure  distribution 
Cp^av.  Figure  15  calculated  from  equation  (59)  may  serve  such  a  purpose 

The  expression  C-jv  =  -2Cp^ay  may  be  considered  as  the  lift  coefficient 

for  the  three-dimensional  wing  tip.  The  rectangular  wing  tip  (k  =  0) 
at  various  Mach  numbers  is  considered.  The  dotted  curves  correspond  to 

(DC  (DC 

—  =  0.2,  and  the  solid  curves,  to  —  =0.5-  The  magnitudes  of  both 

C  ' 

the  real  and  the  imaginary  parts  of  ^  *  decrease  with  increasing 

ae1<m 

Mach  number  as  mentioned  before.  The  change  of  —  from  0.2  to  0-5 

a 

o 

causes  a  slight  drop  in  the  real  part  of  t  .  but  a  greater  drop  in 

otiose 

the  imaginary  part.  The  average  pressure  at  —  =  0.2  deviates  so 

a 

little  from  the  steady  case  (cd  =  0)  that  it  is  hard  to  show  in 
figure  15(a) • 


The  present  method  does  give  a  closed  expression,  but  it  takes  a 
great  deal  of  time  to  calculate  even  a  single  point. 


DISCUSSION  OF  RESULTS 


With  the  above  analysis,  a  few  points  should  be  emphasized.  From 
Evvard's  concept  of  the ' "diaphragm, "  the  two  approximate  methods  have 
been  shown  to  give  a  satisfactory  pressure  distribution  in  the  wing  tip 
influenced  by  the  three-dimensional  flow,  if  the  frequency  parameter  is 
low  and  if  the  leading  and  trailing  edges  are  supersonic. 

The  comparison  of  the  pressure  coefficient  obtained  by  the  second 
method  with  two-dimensional  exact  solutions  of  linearized  theory  is 

coc 

favorable,  if  the  frequency  parameter  — o  is  around  0.5  and 'the  point 
is  at  the  trailing  edge.  It  is  expected' that  the  first  method  should 
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give  satisfactory  results  with  the  frequency  parameter  of  the  order  of 
one-half.  Of  course,  if  only  the  average  pressure  coefficient  over  the 
entire  tip  hounded  by  the  left  tip  Mach  line  is  wanted,  a  higher  fre¬ 
quency  parameter  up  to  the  order  of  unity  may  be  used. 

The  effect  of  the  three-dimensional  flow  on  the  wing  tip  is  essen¬ 
tially  to  reduce  the  pressure  increment  in  the  area.  With  the  approxi¬ 
mation  given  by  the  present  theory,  the  pressure  coefficient  at  the  tip 
side  edge  is  zero  except  at  the  point  on  the  leading  edge.  Of  course, 
the  exact  theory  would  give  finite  values  at  the  side  edge,  which  would 
approach  zero  (as  cd — >0)  as  predicted  by  the  conical-flow  theory.  At 
low  values  of  the  frequency  parameter,  the  effects  of  the  finite  edge 
pressures  give  negligibly  small  contributions  to  the  aerodynamic  behavior 
of  the  wing  as  a  whole.  It  is  in  this  fact  that  the  present  methods  find 
their  justification. 

The  effect  of  increasing  Mach  number  is  to  reduce  the  real  and 
CP 

imaginary  parts  of  v  ,  in  magnitude,  if  the  sweepback  angle  and  the 

oe1^ 

frequency  parameter  remain  the  same.  If  other  conditions  of  the  wing 
remain  the  same,  increase  of  the  sweepback  angle  causes  the  magnitude  of 

CP 

the  real  part  of  “ics’t  increase  and  the  imaginary  part  of  it  to 

ae 

decrease  as  shown  in  figure  12.  Increase  of  frequency  parameter  has 
the  opposite  effect. 

The  pressure  coefficient  of  a  point  influenced  by  the  nose  of  the 
sweptback  wing  is  given  in  terms  of  the  Fresnal  integral.  It  is  not 
difficult  to  carry  out  the  numerical  integration  for  this  case.  In  this 
portion  of  the  wing,  no  approximation  is  used. 

The  first  method  of  calculating  Cp  at  the  tip  using  the  Fresnal 
integral  appears  worthy  of  more  extensive  study  in  the  future  as  the 
method  is  very  simple  and  useful.  The  second  method  gives  closed 
expressions  for  Cp  and  Cp  av  and  ^  ■laas  ^een  investigated  extensively 

in  the  present  report. 

The  present  analysis  applies  only  to  cases  in  which  the  flow  at  the 
leading  and  trailing  edges  is  supersonic,  and  the  nose  Mach  cone  and  the 
tip  Mach  cone  do  not  intersect  within  the  wing  area.  With  slight  modifi¬ 
cation  the  methods  can  be  extended  to  more  general  wing  plan  forms. 

The  Johns  Hopkins  University 

Baltimore,  Md.,  November  15,  1949 
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APPENDIX  A 

SYMBOLS 

a  velocity  of  sound 

A  source  strength 

7  1 

b  =  P7n8  -  —  xT  (equation  (49)) 

7  2 


c(0 


coefficients  of  C 

P,av 

chord 

drag  coefficient 
lift  coefficient 
moment  coefficient 
pressure  coefficient 


Fresnal  integral 


cos  |  (e')2 


f(  )  function  concerning  source  strength 

hjjjg^  maximum  descending  rate 

unit  vector  in  x-direction 


k2' 


constants  expressing  wing  inclination 


6 
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APPENDIX  B 

FORMULAS  USED  IN  SECOND  METHOD 
Formulas  for  Use  with  Equation  (49) 


•2  =  4  -  4im72(x  -  ky)  -  2flD2722(x  -  ky)2  +  ^  im^g^x  -  ky)3  +  icuS^^ 


_  *31  72 
e4  M  7! 


■4im  -  4m272(x  “  ky)  +  2i0337g2(x  _  jjy)  2  +  §  a>S23(x 


= 


6 


732/^2 

M2\71 


-2m2  +  2im ^72(x  -  ky)  +  m^722(x  -  ky)2 


=  !33& Vf2  m3  +  2  ^ 


8=P“fe  3^+Im%(x-ky) 


-  ^34/y2\4  1  -4 
10  "  MM?!  3“ 


%  = 


W^2\l/2 

mH7J 


|  +  |  i^72(x  -  ky)  +  ^  m2722(x  -  ky) 


86  =  mHtI 


-  lif2  ^ 


I  -  I  im72(x  -  ky) 


7  31 

M 


4/72\5/2  _4  1  732 

s8=i?y  “3^ 


x  -  ky)1*- 

-  ky)^l 
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=  2(7i72)ly/2^V2i| 


12 


r31  .+  3  r32 


35 

+  -=^-  700  + 

192  33 


21 

256  73^ 


+ 


1  Zi  +  JL  _L  Zi  -  +  _7_  i  ll  -  +  1  /7l\2~ 

12M2  72  ^  M2  72  192  M2  72  7^2  l60M^\72/ 


*2,0  =  2(7!)l/2(-2) 

ff2,l  =  2(7i)l/2i^72(2  +  |  73i) 

*2,2  =  2(7i)1/25>2722(i  +  I  731  +  |  132  -  ^2 


°2,3  =  2(7i)l/2i^723(-| 


f731-  1^32-  ]§733  + 


I  71  ,  _5_  _1_  Zi-  \ 

6m2  72  48  M2  7 2  73lj 


ct2,4  2(7l)1//2m^72i|'J~i2  "  ¥  731  "  iB  732  "  192  733  "  25^  y3k  + 

+  +JLA.Z lj  +_i_ /Zif 

12M2  72  1+3M2  7 2  31  192  M2  72  32  16oM^\72  /_ 


*3,0 

=  2(71)l/2i^72(73I) 

ct3,1 

=  2(7i)1/2aj2722(731  +  ^ 

:  732) 

a3,2 

=  2(71)1/2ia53723(-|  731 

“S732 

.  35  \ 

288  733/ 

ff3,3 

=  2(7i)1//2^4724(-|  731  - 

2k  732  ' 

.  35  7  _  7  \ 

288  33  128  3k J 

°4,0 

=  2(7l)l/2^2722(|  732) 

CT4,1 

=  2(7)  1///2ia537 23  732  - 

72  ^33) 

CVJ 

-4’ 

=  2(7)1/2a54724(-|  732  - 

I7  . 
72  733 

lfe  734-) 
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=  a(r1)1/2 W  lOi  II  +  _L  U  y  ' 
2’2  1  11  m2  is  r2  W  72  r3iy 


^-3  '  2<U)1/2a“  % 

M 


ili  +  _5.Z± 

12  72  48  7 p 


+  48  ’Z  731  +  —  ~  4 


7  ^1  __ 
192  7 2  732 


Formulae  for  Use  with  Equation  (55) 


a0  "  °l,i 


_  _  io> 

1  ”  U  °1,1  +  2ct1,2 


iCD 


a2  ~  *u~  %2  +  3%3 


„  _  ico  , 

a3  ~  T  %3  +  Wl,4 


±00 


a4  =  —  01,4  +  5ax  5 


ia> 

a5  -  u"  01,5 


b0  =  02,1  +  f%0  +i%0 


b,  =  2o„  „  + 


ico  3 


1  -  ^2,2  tTCT2,1  +  2  03,l  +  °2,1 


b2  3CT2,3  +  T  Q2,2  +  2  ff3,2  +  02,2 


bo  =  4ac  ],  +  a. 


2,4  -  02,3  +  2  a3,3  +  02,3 


,  ico 
°4  ~  u  02,4 


„  _  _  '  ia>  5 

c0  "  03,1  +  ~  a3,0  +  2  \o 


1  _1_ 
180  m2 
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io)  7 

x0  -  °9,1  +  u  a9, 0  +  2  aio,o 


10) 

xi  "  u  a9,i 


J  =  cr 


10,0  u 


Formulas  for  Use  with  Equation  (55)j  a)  =  ^ 

p  a32 


1,1  "  3 


°i,s  =  jH-21) 


_  1  =2/  2  1 
°1,3  “  3  m  '3  ‘  ,u2 


"P  6  T  4m2 


i  im3fi  +  — 


1  =4/1  ^  l  ,  l 

=  7T  00  hrs-  +  - rr  +  - 


al,5  3  “  \30  10M2  8oi? 

CTn  n  =  3-1/2(-4) 


o21  =  3"1/2  iai(2) 


a2  2  =  3-1/2  ai2(f  - 
d>d  \3  3M2 


ap  o  =  3"1/2  ia)3  (-  i  + 

2,3  l  6  4M2 


_  0-1/2  -4/  1  ,  1  L  1 

2,hr  m  (  30  10M2  801?, 


a3,0  =  P'l/2  m-J 
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_  R-1/2  =2/  8 

03^1  =  P  '  CD  ( -  — 


-1/2  iffitfi 


“3,2 

a^-p-VSstfi 

CT4,o  =  P"l/2 


i+5 


ff4,l  =  P' 


•V2 


=  P"l/2 


ct5,0  ~  p 


105/ 


5,1 


=  p-1/2  sV  JU 


V525/ 


%o  -  ^  *W 


„  ni/2  1  =2/2 

a7,l  =  p/  P“(l 


Or,  ,  =  61/2  J,  3=3  ^l 


7.2 

7.3 


M‘ 


a„  ,  -  tt1/2  i  aV.i  . 

M‘ 


2  “  V5  ' 


-  r1/2  1  =2/8 

°s,o  -  p  ^“(9, 


“8,1  -  e1/£  ^  ^(-f 


%*  =  »1/2  hM-h  -  - 


30M' 
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(M>  I) 


5* —  Pressure  di stribution  along  trailing  edge  of  a  rectangular 

wing  of  unit  chord  with  ^  =  0.2,  M  =  ^2,  and  x  =  1.  Calculated 
with  first  method  (Fresnal  integral). 


(M>l) 
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;ure  6.-  Wing  area  Swl  (PABC)  which  affects  the  potential  at  P(x,y) 

and  the  oblique  coordinate  system  (x*,y* )  chosen  for  calculation  in 
second  method. 


(a)  Real  part  of 


Figure  7*-  Real  and  imaginary  parts  of 


along  -y  at  various 


values  of  x  for  a  rectangular  wing  tip  (k  =  0)  of  unit  chord  with 
—  =  0.2  and  M  =  \J2  (used  for  calculating  pressure  distribution 
on  wing  tip  influenced  by  three-dimensional  flow).  Second  method. 


-y 


(to)  Imaginary  part  of 


rtCp 

aeicut' 


Figure  7--  Concluded. 
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(a) 


-y 


(a)  Real  part  of 


nC- 


ae 


last 


(b)  Imaginary  part  of 


itCn 


ae 


icat 


jcC 

Figure  9.-  Real  and  imaginary  parts  of  along  -y  at  various 

aeimt 

values  of  x  for  a  rectangular  wing  tip  (k  =  0)  of  unit  chord  with 

CDC 

—  =  0.2  and  M  =  3. 
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(a)  Real  part  of 


ae 


ia)t 


(b)  Imaginary  part  of 


ae 


icot* 


Figure  11.-  Comparison  of  results  of  second  approximation  with  two- 

dimensional  solutions  at  left  tip  Mach  line  of  right  rectangular  wing 

tip  (k  =  0)  along  chord  at  =  0.283  and  =  0.707  with 

a|32  a|32 
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(a)  Real  part  of  - E— . 

a.eicuk 

(b)  Imaginary  part  of  — 

ae 


ia>t' 
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Figure  12.-  Comparisons  of  results  of  second  approximation  with  two- 
dimensional  solutions  at  left  tip  Mach  line  of  a  right  sweptback  wing 
tip  along  chord  for  various  values  of  k.  k  =  tain  A  =  0.25  anil  0.  5. 

ODC  ^  ^  ' 

— -  =  O.283,  and  M  =  t/2. 
a3 
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(a)  Real  part  of 


ae 


lost 


C 

Figure  13.-  Distribution  of  — -E-  along  constant  x  -  y  at  various 
D  10Jt 


ae 


CDC  i 

values  of  x  +  |By  for  a  45°  sweptback  wing  tip  at  M  =  2,  —  =  0.4, 

a 


10 
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Figure  13.-  Concluded. 
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Figure  l4. -  Comparison  of  results  of  second  approximation  with  two 
dimensional  solutions  at  left  tip  Mach  line  of  a  sweptback  wing 

tip  along  chord.  Sweepback  angle.  4-5°;  M  =  2;  —  =  0.4.  k  =  1. 


(a)  Real  part  of 


aeic0t* 


Figure  15- -  Average  pressure  distribution  over  entire  rectangular  wing 
tip  (k  =  0)  influenced  by  three-dimensional  flow  as  it  varies  with 

Mach  number  for  ^  =  0,  ^  =  0.2,  and  ^  =  0.5. 


BSaSM 


numbers,  sweepback  angles,  and  frequency  parame-  numbers,  sweepback  angles,  and  frequency  parame¬ 
ters.  Comparisons  of  pressure  coefficients  calcu-  ters.  Comparisons  of  pressure  coefficients  calcu¬ 
lated  from  present  methods  and  exact  two  dimen-  lated  from  present  methods  and  exact  two  dimen¬ 
sional  solutions  based  on  linearized  theory  are  made.  sional  solutions  based  on  linearized  theory  are  made. 
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